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ABSTRACT 
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Detailed analysis of short-period travel time, dT/dA and 
waveform data reveals the upper mantle structure beneath an oce- 
anic ridge to depths of 900 1cm. Mote than 1400 digital seismo- 
grams from earthquakes in Mexico and central America recorded 
at SCARLET yield 1753 travel times and 5B direct measurements of 
dT/d A as well as high-quality, stable waveforms. The 29 events 
combine to form a continuous record section from 9* to 40* with 
an average station spacing of less than 5 km. First the travel times 
are inverted using the lau method of Bessonove et al. (1974,1976); 
the resultant model is adjusted to agree with the experimental p-A 
values. Further constraints arise from the observed relative 
amplitudes of mantle phases, which are modeled by trial-and-eiror 
using WKBJ synthetic seismograms (Chapman. 1976; Wiggins. 1976). 
Model GCA. which is appropriate for western Mexico north of 20* 
latitude, is similar to existing upper mantle models for shield, 
tectonic-continental, and arc-trench »-egimes below 400 km. but 
differs significantly in the upper 350 km. CCA velocities are very 
low in this depth range; the model "catches up" to the others with 
a very large velocity gradient from 225 km to 390 km. This well- 
resolved feature is consistent with the sheai^wave model TNA for 
western North America of Grand k Helmberger (1983). The abun- 
dant data from 20* to 30* control the detailed shape of the 860-km 
discontinuity. Very high velocity gradients lie both above (620-880 
km) and below (861-680 km) a 2.8% velocity change. 


Introduction 


The nature of lateral variations In upper mantle seismic velocities is a prob- 
lem of broad geophysical interest. Velocity structure determinations for many 
regions can place constraints on both the seals of mantle convection and the 
depth extent of velocity differences between continents and ocean basins. 
Hager k Raefsky (1961) predict large depressions of a chemical “8iU km" 
discontinuity beneath subducted slabs if convection is confined to the upper 
mantle. Sipkin k Jordan (1975, 1978) suggest that lateral differences between 
shields and old Deans extend to 400 km depth to satisfy multiple ScS travel 
times, while Okal k Anderson (1975) insist that most of the differences between 
shields and old oceanic ScS data are explained by heterogeneity shallower than 
200 km. Recei ’ upper mantle models for continental shields (e.g. King k Calcag- 
nile, 1976; Given k Helmberger, 1980), "young" continental regions (Johnson, 
1967; England et al.. 1977; Burdick k Helmberger, 1978; among others) and 
island arc regimes (Kanamori, 1967; Fukao, 1977) have been constructed using 
body waves; these models tend to converge below 200 km. Very little work has 
been done, however, on the important areas of ocean basins, continental rifts 
and spreading centers. Some surface wave studies indicate low upper mantle 
velocities for young oceans (Knopoff et al.. 1970; Montagner k Jobert. 1981; 
Wielandt k Knopoff, 1582) to depths of at least 200 km. England ct al. (1978) 
analyzed P-wave data from the North Atlantic Ocean, and Green (1978), Nuiet k 
Mueller (1982) and Lenartowicz k Albert (i960) studied the African rift region, 
for which very slow teleseismic travel times have also been documented. The 
detailed characteristics of the upper mantle to 1000 km beneath spreading 
centers, however, are unknown. 

We have investigated the upper mantle P-wave velocities under the Gulf of 
California spreading center. Mexican earthquakes recorded at the California 
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Institute of Technology - U.S. Geological Survey Southern California Seismic Net- 
work (SCARLET) provide a unique, dense, high-quality data set. The narrow 
azimuthal range of the epicentral distribution results in a nearly ideal data 
profile. We exploit the large amount of travel-time, apparent velocity (dT/dA) 
and wa'-eform data in the modeling process. While the travel times control the 
model's gross integral properties, the dT/dA measurements provide information 
about the absolute velocities at the rays' turning points. The relative ampli- 
tudes of phases are most sensitive to the velocity gradients near the bottoming 
points. We combine these data by first inverting the travel times, perturbing 
that model to fit the p - A data, and then performing trial-and-error synthetic 
seismogram modeling to fit the short-period waveforms. The final model 
satisfies all three data sets. Thus our model has more resolvable detail than 
those derived from travel timei. alone (e.g. Hales, 1972; Mass^, 1973, 1974; 
Green. 1978) or from array studies utilizing travel times and apparent phase 
velocities but which nntain no waveform analysis (Johnson, 1967; Simpson et al., 
1974; Ram k Mereu, 1977; King k Calcagnile, 1976; Ram et al., 1978; England et 
al., 1977, 1978; and others). In addition, our high-qual ; ty, dense array data 
allows better resolution than discrete source-receiver waveform studies such as 
Helmberger k Wiggins(l971), Wiggins k Helmberger (1973), Dey-Sarkar k Wiggins 
(1978) and McMechan (1979). While short-period waveforms are not as stable as 
the equivalent long-period data, teleseismic waveforms recorded across the 5° 
aperture of SCARLET are very reproducible, indicating good stability for simple 
events. Our final model represents a synthesis of differing constraints and 
results in o well-resolved, detailed view of the upper mantle under or near an 
active oceanic ridge. 
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The Data Set 

The scismically active area6 of the Gulf of California, Rivera Fracture Zone, 
East Pacific Rise, and Middle America Trench are the source regions for this 
study. Figure 1 illustrates the experimental geometry, including 2* arcs drawn 
at the travel path midpoints of the 22 events at distances of less than 30*. 
Clearly, the upper mantle sampled by these earthquakes is not influenced by the 
Middle America Trench, but represents the Gulf of California and adjacent exten- 
sional areas. The events occurring on the spreading centers and fracture zones 
have strike-slip mechanisms oriented unfavorably for P-wave radiation to SCAR- 
LET. and tend to have complicated source signatures. Most of these events are 
closer than 20*. The sutduction zone earthquakes, on the other hand, at epi- 
central distances greater than 18* , are dip-slip events, which produce ample P- 
wave energy and are often simple in character. The events range In distance 
from 9° to 40° and occur in the narrow event-station azimuth band of CIO* to 
345°. Varying in depth from 10 km to 150 km, they have body-wave magnitudes 
of 5.0 to 6.3. All events occurred between September, 1977 and December, 1979, 
and are listed with the PDE epicentral information in Table 1. 

Each earthquake is recorded by the short-period vertical, digital, triggered 
CEDAR s/stem (Johnson, 1979) at the California Institute of Technology. In its 
current configuration SCARLET has more than 200 stations; from 1977 to 1979 a 
well-recorded teleseism would trigger 120 stations, about 60 of which fit the cri- 
teria for inclusion in the data set (Figure 2). Elongate in the northwest- 
southeast direction, the array has an aperture of 5* and irregular station spac- 
ing averaging 25 km. Although the array stations have varying instrumentation, 
the responses are very similar at 1 Hz, the predominant frequency of the tclese- 
ismic signal. 
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Over 1400 digital seismograms were collected from the 29 events, yielding 
1753 travel times Including 438 which are secondary arrivals. In addition, we 
obtained 58 direct measurements of the ray parameter, d?/dA, spanning the 31° 
distance range. Figures 3, 4, 5. and 8 show examples of record sections for 
several events at different distances. Because SCARLET is not well-calibrated, 
only relative amplitudes are used, and each trace is scaled to its maximum 
amplitude. Adjacent records are very similar, indicating excellent waveform 
stability across the array. The events in Figures 4, 5. and 6 are all simple and 
impulsive, allowing unambiguous selection of secondary phases. Each record 
section covers 4* - 5* in distance and collapses about 10* of azimuthal variation 
onto a plane. Distances in Figures 3-6 are not corrected for event depth. Por- 
tions of interesting upper mantle triplication phases are visible for each event, 
but a more complete picture is obtained by combining the 10 cleanest events 
covering the entire distance range into oi.e record section (Figure 7). This 
representation contains 477 depth-corrected seismograms with an average data 
spacing of B km. Inclusion of all available data reduces i.he spacing to less than 
5 km. Secondary arrivals from both the "400 km” and "670 km" discontinuities 
are seen clearly from 14* to 28“ . This high-quality data prompts careful and 
complete data analysis to insure a robust, detailed upper mantle model. 
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llala Preparation and Analysis 

Receiver structure in southern California 

An area of complex geology and present-day tectonic activity, southern Cal- 
ifornia has a complicated, heterogeneous crust end uppermost mantle structure 
which affects incoming teleseismic signals. The lar^e amount of available data 
(both local and teleseismic l„) have prompted several studies probing the 
nature of the receiver structure beneath SCARLET. Kanamori k Hadley (1975) 
reported on the region's upper crustal velocities; Lamanuzzi (1981) and Hearn 
(1983) have investigated gross crustal and upper mantle properties using Pn 
travel times. Teleseismic P arrivals were used in upper mantle heterogeneity 
studies by Hadley k Kanamori (1977), Raikes k Hadley (1979), Raikes (1980) and 
Walck A Minster (1982). The spatial pattern of the strong azimuthai variation of 
teleseismic P-resi duals (Raikes, 1980) is consistent with a high-velocity body in 
the upper mantle beneath the Transverse Ranges in southern California, first 
proposed by Hadley k Kanamori (1977) and verified by W'alek k Munster (1982). 

Correcting for near-receiver velocity variations in a multi-azimuth data set 
in this complex area could require detailed ray tracing, but for our "single" 
azimuth data e simpler approach was adopted. We constructed station correc- 
tions derived from travel times of more distant (30* < A < 40*) central American 
earthquakes (Figure 1. Table l) in the same azimuth band. These events are 
free of complicating upper mantle phases, yet the rays are incident et the 
receiver at angles similar to the closer events. Planes are fit to arrival times 
using least-squares for several large, impulsive events; the station coi rections 
are the averaged station residuals from the plane-predicted arrival times. The 
procedure assumes that 1) constant corrections are appropriate over the entire 
30* azimuth spread and 2) the travel-time curve is smooth beyond 30* in 
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distance. 

The seven distant events have a total azimuth range of only B° and are 
remarkably consistent: 96 stations have average residual values with standard 
deviations or less than 0.10 s (shown in Figure 2). Our empirical corrections 
include effects of both local structure and elevation, and are applied to all the 
data to reduce travel-time scatter and aid in identification of secondary phases. 
An example of a record section before and a er application of the emp'rical 
corrections is shown in Figure B. 

w e might expect a close correlation between the empirical adjustments and 
Raikes' (1980) teleseismic residuals for the same azimuth range. Although her 
data are from more distent (A ~ 55*) earthquakes- and are single station residu- 
als instead of deviations to the array least-squares plane, the contour plots of 
the two residual sets (Figure 9) are similar in shepe. The difference in absolute 
magnitude of the residuals occurs because Raikes* (1980) residuals are refer- 
enced to an individual station, CSC, instead of the plane average and have been 
corrected for elevation, basin sediments and large-scale Moho depth variations. 

Travel Times 

The 29 events are retrieved from magnetic tape storage and the travel 
times hand-picked with accuracy that approaches the digitization interval. .02 s. 
Both direct picking and cross-correlation techniques were tested with nearly 
identkui results. The results displayed here are for hand-picked times. Many 
records are low-gain or noisy; about one-half of the original seismograms are 
ultimately rejected, leaving about 60 records for each event. Each travel time is 
corrected for ellipticity, for depth (using the Jeffreys upper mantle model) and 
with the empirical station correction. The 1753 travel times provide a continu- 
ous curve from 9* - 40° (Figure 10). Errors in the earthnuakes* hypocenters and 
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origin times still cause considerable scatter in the travel time data. These 
uncertainties ere removed by applying baseline shifts based on the average of JB 
residuals in completely overlapping 0.5* distance windows. These shifts preserve 
the shape of the observed travel time curve, yet reduce the data scatter to 0.4 s; 
the agreement In differential travel times f^r secondary phases is excellent (Fig- 
ure 10b). 

Apparent Velocities 

The numerous high-quality travel times allow calculation of many reliable 
dT/dA, or ray parameter, estimates for both primary and later arrivals. The 
standard plane-fitting technique (e.g. Manchee k Weichert, 1968) for phase velo- 
city determination is used on the empirically corrected but not baseline-shifted 
times. Because of the great size of the array, the earth’s sphericity is taken 
into account according to the procedure outlined by Walck k Minster (1982). 
Before computation begins, each event record section is visually inspected for 
changes in the travel-time slope with distance. When such a change exists, the 
array is divided into two sections and two ray parameter estimates are made. 
Measurements of dT/dA made with only part of the array are still stable and 
accurate because of the network’s large number of stations. 

We obtained 58 p - A points (Table l) shown in Figure 11. Fourteen of the 
measurements are determined from later arrivals. The values have low e’Tor 
estimates, are stable and exhibit little scatter. Some smoothing of changes in 
apparent phase velocity could occur over SCARLET'S 5* aperture. To test *or 
this, the network is split in half by a rortheast-aouthwest trending line and all p - 
A points redetermined and compared to the whole-array estimates. We found no 
significant differences between the 77-point split-array data set and the 58-pomt 
whole -array group. The excellent travel-time and dT/dA data permit standard 
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inversions to find a reasonable sorting model for synthetic seismogram model- 
ing of the structure. 

Inventions 

Travel times 

Instead of inverting the p - A curve with the classic Wi iert-Herglotz for- 
mula, we performed a direct travel time inversion using the tau method of Bes- 
sonova et al. (1974,1976) which has advantageous statistical properties. The 
quantity 

r(p) = T-pX 

where T is travel time, p is ray parameter and X is epicentral distance is calcu- 
lated directly from the travel time data. For a fixed ray parameter, p,, r(p 0 ) Is 
the extremum of the function T(X) f=Pt along a single travel tune branch (Besso- 
nova et al.. 1974, 1976). To statistically determine t(j> u ). we assume that 
T(X)p mPt is a constant for a small window in X centered on X 0 , the distance at 
which t(X) is an extremum. Then r[p 0 ) is a simple average ot all the r(X) pmpt 
data points in that X interval (see Figure 5 in Bessonova et al., 1976). A 
confidence interval, <J # , for the estimate is given by 

6, = n~ l/9 s t m (n— l) 

Here n is the number of observations, s is the standard deviation cf r(p 0 ) and 
t m (n-1) are values for the Student's distribution with n-1 degrees of freedom 
(see Bessonova et al., 1976 for details). The statistical confidence limit is impor- 
tant in assessing the model uncertainty. 




\ 
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The Mexico data set yields 33 values for r in the p range 0.3 s deg to 14.0 
s deg _l , which roughly corresponds to distances of 10* to 40*. Typical values 
for n and are 50 and .20 s, respectively, although some values are much 
better determined than others. After assuming a crustal model, the tau data 
was inverted first for the best single velocity model, and then, incorporating the 
uncertainty values, for the extremal bounds allowed by the travel times (Figure 
12). The tau method does not utilize our independent p - A data and thus it may 
produce a model that fits the travel times but not the experimental dT/dA 
curve. By perturbing the individual r(p) values within their statistical uncer- 
tainties, we obtain a model whit, matches both the travel time and apparent 
phase velocity data. The mode, \ , gure 12 has no low-velocity zone because of 
a lack of compelling evidence for one in the waveform data. A 32 km thick crust 
was used arbitrarily; a 15 km thick oceanic crust may also be appropriate to 
represent the model area. The somewhat gradual crust-mantle transition is 
required by the very slow observed travel times from 9* io 13° (see Figure 10). 
The gradient above the 400 km transition zone is controlled by one group of 
lower-quality data points which will be discussed in more detail in a later sec- 
tion. The shape of the extremal bound envelope indicates that the bost- 
determined model depth ranges are from 450 - 825 km and below 700 km. This 
tau model is the starting point for the Lrial-and-i rror waveform inversion using 
synthetic seismograms. 

Waveforms 

To make a synthetic seismogram, commonly a series of convolutions are 
performed in the time domain: 


y(t)=s(t)*m(t)*a(t)*i{t) 
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where y is the seismogram, s ir the source-time function, a is the attenuation 
operator, i is the instrument response, and m is the Green's function for the 
travel path (Helmberger & Burdick, 1979). In order to isolate the desired quan- 
tity, m(f), estimates ofs(f), a(f). andi(f) must be made. 

The source-time functions *(t) for events of m* ** 5.5 are often short. (2-3 
s) and impulsive, especially for subduction tone earthquakes. Thus for many of 
tue events, a clear pulse winch bottomed in a smooth portion of the mantle is 
the first arrival, separated from and followed by the reflected or refracted 
mantle-generated phase (see Figure 4 for an example). In this case the fl r st 
arrival -epresents the source-lime function s(t) convolved mth the instrument 
response i(t); this is valid acrors the entire array. The attenuation operator is 
neglected; possible problems associated with this are discussed belor. By 
extracting an high-quality flrst arrival from an tual digital seismogram to use 
us the source-instrument response, we eliminate source uncertainty frorr the 
modeling problem. Cf course each event must be modeled with the appropriate 
source wavelet. Only one convolution is made in the synthetic calculation: man- 
tle response with the source vravelet time ser as. Only events with simple 
source-time function* are U3ed in the modeling process. 

To compute the nantle Green’s functions, the WKBJ method (Chapman, 
1976; Wiggins, 1978) is utilized. Formulated for inhomogeneous media, this tech- 
nique cannct readily accommodate causal attenuation and is not as accurate 
near first-order discontinuities and for grazing incidence as some other 
methods. It is. however, very rapid to compute WKBJ seismograms so that a 
wide suite of models can be tried at little expense. This modeling allows the 
inclusion of relative amplitude and waveform data which are sensitive to velocity 
gradients in the mantle. It is used to "fine-tune" whe models already derived 
from travel time and dT/dA data. 


We cbose a subset of seven simple events to model, covering the distances 


13* to 30°. Starting from tne tau model derived earlier, relative ampiit'ide (aid 
timing problems were identified and corrected. Then the T * A and p - A plots are 
generated for the revised model to insure that those parameters remain accept- 
able. This iterative scheme is continued until all three data types are matched. 
The model is then simplified as much as is possible while still matching the data. 
The final model, GCA, is the result of more than 100 iterations using waveforms. 
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Results 


Relative Amplitude Pattern* 

Model CCA (Figure 13). which inc orpomtes vavr.fonn data, differs consider- 
ably in detail from the trawl-time derived starting model (Figure 12). Relative 
amplitudes provide constraints which eliminate many models allowed by the 
travel time and imperfect ray parameter measurements. In our data set, many 
separate events overlap in distance and are characterized by similar and stable 
relative amplitude patterns, which we use in the synthetic seismogram model- 
ing. 

. From 9* to 13* , event signati’-es are generally complicated and exhibit less 
consistency array-wide than do the more distant events. The two events in this 
distance range (nos. 3 and 10 in Table l) are small, with m*s of L.O and 5.3 
respectively. As they are located in the Gulf of California, their mechanisms are 
probably strike-slip with a P-wave node facing the array; the source characteris- 
tics ere complex, va r y with azimuth and so often result in poor records. While 
first arrivals are small, no weakening trend with increasing distance or obvious 
secondary arrivals argue conclusively for the presence of a well-defined low- 
velocity zone. Travel times are very slow in this distance * ange, and apparent 
velocity measurements are scattered and less reliable due to poo^ signal-to- 
noise ratios fer many records. 

Very weak initial arrivals and energetic secondary phases characterize 
seismograms at distances of 14* to IB* (see Figure 3). At 15*. the secondary 
arrival, which is the reflection fro /i the 390 km discontinuity, is 6 s behind the 
first break. Complex interference patterns develop near 1B # as the later phase 
moves through and takes over as the first arrival. 
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Small near 19* , the first arriving energy increases in amplitude relative to a 
strong secondary phase until a distance of 21 # . Now the initial-arriving energy is 
bottoming in the smooth portion of the mantle between 400 and 650 km. At this 
range both wavelets are equal in amplitude; the second wavegroup arrives 3 to 4 
s behind the first (Figures 4, 5). The first arrival weakens near 23* and an 
interference pattern emerges at 24° between the two phases. This amplitude 
pattern is duplicated for severed events in this distance range. The strong later 
arrival observed here is the part of the EF branch (Figure 10b) of the travel 
time curve which bottoms at the lower discontinuity. For many u* per mantle 
models (e.g. those of King k Calcagnile, 1976; and Burdick k Helmberger, 1978) . 
the AB branch would be moving out in this distance range with large amplitudes 
and would be highly visible. The null-observation of this travel time branch 
plays an important role in the modeling process. 

From 26° to 28° a secondary arrival moves out and weakens relative to the 
first arrival (see Figure 6). This is the "back branch" rf the second triplication 
and u part of the CD travel time branch. Seismograms become simple in char- 
acter at distances beyond 28°. 

Model Description 

The crust, uppermost mantle and low-velocity zone of GCA (Figure 13. Table 
2), while subject to some bounding information, are non-unique. Above 125 km. 
the model parameters are adjusted mainly to agree with travel times, subject to 
several constraints: a single layer, 20 km crust is a compromise between the 
30-km, two-layer southern California structure and a thin oceanic crust. The Pn 
velocity is set at 7.9 km/s, the best value for southern California (Hearn, 1983) 
and assumed to be appropriate for continental areas of \orthwest Mexico 
covered by GCA. Travel times beyond 13° control the integral of the size of the 
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low-velocity rone and the absolute velocities above 125 km. 

Arrivals closer than 13.5* are in the shadow zone of GCA (see Figure 10b). 
An offset in the travel time data is also observed at that distance. Fig ure 1 
■hows that the two closest events, in the Gulf of California, have midpoints 
beneath the Gulf itself, while the events from 13* - 16* (nos. 2.4,12.13 in Table l) 
are located on the Rivera Fracture Zone, and their rays turn under the Baja Cali- 
fornia peninsula. The offset travel times could be due to strong lateral varia- 
tions between the two regions for depths of less than 150 km. Travel times from 
only events 3 and 10 (Table l) and an arbitrary crustal structure were used in a 
separate tau inversion for shallow structure of the Gulf itself. Figure 14 shows 
that, in the absence of a low velocity zone that would put 9° - 13° in a shadow 
zone, the transition from crust to mantle must be gradual for the spreading 
center (Figure 14. Table 3). A smooth, gradual crust-mantle transition is also 
documented for the Jordan-Doad Sea rift (Ginzburg et al.. 1979). The travel 
time? oT Figure 10b, then, can be satisfied by model GCA' (Table 3) for distances 
or less than 13° and GCA beyond 13.5*. While two models are proposed based on 
the shift in travel times at 13*. both GCA and GCA' satisfy the waveform data 
from 13* to 15* reasonably well. These two models are based on the differing 
crust-mantle transitions between the Gulf of California and adjacent continental 
areas. While the Gulf itself may have no seismic lid. a gradual crust-mantle tran- 
sition and no velocity reversal, the continental portion of the study area, with 
more usual Pn velocities, requires a region of negative velocity gradient to 
satisfy travel time data. 

Tight constraints on the model shape begin at depths of 125 km. The small 
amplitude first arrivals from 15* to 17* require a very slight positive velocity 
gradient between 125 and 225 km (see Figure i5). A model for the western 
United States. T7 of Burdick A Helmberger (1978), has a more moderate 
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gradient in this distance range and produces large first arrivals at i5* (Figure 
15). The first arrivals of GCA are still large relative to the 390-km reflection at 
these distances, but Q may have an important effect here. The initial arrivals 
spend more lime in the highly-attenuating asthenosphere than du the mantle 
reflections. 

A flrst-order discontinuity of 4.9* at 390 km produces large amplitude 
secondary arrivals at 14* which become first arrivals near 18° (see the CD 
branch of Figure 10b). All discontinuities in GCA are represented as steps in 
velocity because equivalent gradients over 10 - 20 km are not resolvable. In 
many regions, the "back branch" of the 390 km travel-time triplication (AB 
branch) is observed to distances of 24* (England et al., 1977; Burdick A Helm- 
berger, 1970; and others) or even past 30* (King 6c Calcagnile, 1976). This is 
indicative of a small velocity gradient between 300 and 400 km which is 
inefficient at turning energy to the surface so that it is seen at larger ranges. 
The data used in this study show no evidence for the AB branch past 20° for 9 of 
10 events in that distance interval. The anomalous event has a complicated 
source and a low signal-to-noise ratio, and occurred on a fracture zone rather 
than in the subduction regime (Figure 16). In the tau inversion, dubious secon- 
dary times from that event (no. 28 in Table l) are used to help define the first 
discontinuity, so the resulting model has a more modest gradient above 400 km. 
Synthetic seismograms for this structure, however, fit observations from 20* to 
23* very poorly (see Figure 17). The absence of an observed AB branch past 20* 
thus requires the steep velocity gradient seen in GCA from 225 to 390 km. This 
unusual gradient is well-supported by the data and is a feature significantly 
different from models proposed for shields, trenches or tectonic-continental 


au’eas. 
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Large first arrivals are produced from 20° to 23° by the strong velocity gra- 
dient from 390 fo 620 km depth; a small inflection at 540 km amplifies initial 
pulses near 21* (Figure 17). The shape of the velocity-depth curve from 820 to 
700 km is particularly well-resolved by large amounts of data recorded in the 
appropriate distance interval, 22* - 28*. Figure 18 shows synthetic seismograms 
generated for different shapes of the 680 km discontinuity compared to the 
data. The very fast velocity increase from 620 to 660 km is necessary to gen- 
erate the correct relative amplitudes between the EF branch (first arrival) and 
CD branch (second phase} on the seismograms at 25.9* and 27.1*. Improved 
amplitude and timing relationships at 23 # are gained by increasing the gradient 
just below the discontinuity. Thus a first order velocity jump of only 2.8%. cou- 
pled with large gradients immediately above and below, satisfies the waveform 
data best. 

Observed seismograms are simple beyond 28" in distance (Figure 6). The 
gradient which fits the p -A data is adopted in CCA and is very similar to that for 
a JB earth. 

The need for a sharp velocity gradient from 250 to 390 km and the 
increased resolution of the fine structure of the 660 km velocity break mark the 
improvements in detail of the waveform-constrained model GCA over the earlier 
tau-derived model for the same data set. The addition of waveform techniques 
to traditional array analysis of short-period data are important in elucidating 
the fine structure of the upper mantle. 
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Discussion 

Many past studies of upper mantle structure have relied on less than high- 
quality short-period travel time and waveform data, t ; tber from discrete 
receivers or small arrays. The advantages of large-aperture arrays such as 
SCARLET or NORSAR are obvious: stable dT/dA measurements can be made 
using the whole network or subsets of it, secondary arrivals with differing phase 
velocities are easily identifiable, and anomalous traces, such as those contam- 
inated by unusual receiver structure, can be identified and discarded. For suit- 
ably simple events, waveform modeling is a usef-J tool for short-period data as 
well as the more stable long-period energy. As more high-quality digital data 
becomes available, additional detailed investigations including many data types 
should result in less ambiguous, better resolved upper mantle models. 

A comparison of the spreading center model GCA with well-constrained 
models for differing tectonic regimes is very revealing. We consider a continen- 
tal shield model KB (Given A Helmberger. 1980), tectonic-continental model T7 
(Burdick Sc Helmberger, 1978), and trench-arc model ARC-TR (Fukao, 1977) (Fig- 
ure 19). K8. which represents northwest Eurasia, was derived from synthetic 
seismogram modeling of both long- and short-period P waves and is constrained 
to fit the NORSAR-determined p - A curve of model KCA (King & Calcagnile, 1976). 
Similarly. Burdick Sc Helmberger'b (1978) western United States model T7 relies 
on the Johnson (1967) apparent velocity measurements (made at TFSO in 
Arizona) as well as waveform modeling of earthquakes with known source 
mechanisms. An unusually complete travel time and dT/dA data set constrain 
the parameters of ARC-TR, a model for the Pacific Ocean trench near Japan. Fig- 
ure 19 shows that all four velocity-depth curves are quite similar below 400 km 
depth, with some slight differences in velocity gradient at depths greater than 
670 km. The small disparities in the depths to the "400 km" and "670 km" 
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discontinuities are probably not resolvable within Ihe data constraints and 
modeling error. 

For depths shallower than 400 km, however, significant differences between 
the regions appear. In all the models. Lid and low-velocity zone shapes are non- 
unique, with the arc and shield models predicting considerably faster arrival 
times for regional distances than do the young continent and ridge representa- 
tions. At about 200 km, K8, T7 and ARC-TK converge, while GCA features much 
lower velocities. Low (ARC-TR) to moderate (KB. T7) velocity gradients prevail 
between 200 km and 400 km except for GCA. in which velocity increases very 
rapidly with depth in that depth range; GCA velocities mrrge with the other 
models at 350 kin. 

Low velocities for both P and S waves are well-documented for oceanic 
ridges and continental rills. Surface wave dispersion studies require a 2 ry low 
S-wave speeds beneath young ocean (e.g. Knopoff et ad.. 1970; Montagner k 
Jobert, 1981; Wielandt 6c Knopoff, 1982). Detailed refraction work performed on 
rifts (Ginzburg et aL, 1979; Maguire A Long, 1976; Puzyrev et aL, 1973) supports 
low values for Pn velocities. Oceanic ridge refraction studies (e.g. Gettrust et al. 
1982; Lewis 6c Germany, 1982) find compressional speeds of 8.0 km/s at shallow 
(3 km) levels, but this "lid" may be very thin (see Bulin, 1979) and underlain by 
extremely slow material. 

There is little data pertinent to the deep structure of rifts and ridges. Avail- 
able P-residuals (Rowlett k Forsyth, 1979) and PP-residuals (Dorbath k Dorbath, 
1981) for the Mid-Atlantic Ridge are large and positive, indicating very low veloci- 
ties, possibly to great depth. Very late P-wave arrivals at Addis-Ababa (Dziewon- 
ski k Anderson, 1983) are also observed for the East African Rift. Thus the velo- 
city value of only slightly over 8 km/s at 200 km in GCA is consistent with the 


available data. 
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Smali observed first arrival amplitudes from 14* - 16° require a very slowly 
increasing velocity from 100 to 200 km. The gradient, however, trades off with a 
possible low Q zone just below the lid which couli reduce the amplitudes of the 
first arrivals near 14°. While such an attenuating zone will probably have a 
greater effect on wave amplitudes than on periods, the synthetic seismograms 
indicate that there is no noticeable frequency depletion of the first arrival rela- 
tive to the second at 14*. 

The main reason for GCA’s steep gradient from 225 - 390 km is that no 
arrivals corresponding to the AB travel time branch (Figure 10b) are observed 
beyond a (surface focus) distance of 20°. In some shield regions (King & Calcag- 
nile, 1976), this branch is noted past 30°. Typical observational limits for tec- 
tonic continental areas are 24° to 26° (Wiggins Sc Helmberger, 1973; Ram et al., 
1977; Burdick Sc Helmberger, 1978). England et al. (1978) used oceanic events 
recorded at NORSAR for their model NAT and see the AB branch to 24 , 1 but for 
that range the rays’ bottoming points are no longer beneath young ocean, so a 
different structure migut be expected. In Australia, Simpson et aL (1974) do not 
observe an AB branch beyond 21*. Model SMAK 1 has a small velocity gradient 
above 400 km which predicts AB arrivals well past 30°; Simpson et al. appeal to 
a properly placed low Q zone to suppress amplitudes of the AB phase. For Indian 
Ocean earthquakes recorded at the Gauribidanur array in southern India, Ram Sc 
Mereu (1977) cannot identify the AB branch past 19*. Their model RM-J also has 
a shallow gradient above 400 km but terminates the AB branch with a very deep 
(175 - 332 km) low velocity zone. 

Other studies with oceanic sources, therefore, document the absence of the 
AB branch beyond 20* but resort to non-delerministic methods to diminish that 
phase's amplitudes. Since low upper mantle velocities are consistent for »‘ ; dges, 
the very slow uppermost mantle of GCA underlain by a region of unusually fast 



velocity increase with depth is very appropriate and a less artificial way of 
effectively shortening the predicted AB travel time branch. 

Shear wave data from the East Pacific Rise support the idea of a large velo- 
city gradient between the depths of 2C ^ and 400 km. Grand & Helmberger's 
(1983) model TNA (Figure 20), derived from long-period SS phases recorded at 
North American stations, is very similar to GCA in general character, and 
includes a high gradient from 250 - 400 km in depth. 

The discontinuities near 400 km in the four models of Figure 19 are all quite 
consistent in size and shape. The absolute depth to the velocity transition is 
dependent on the assumed shallow structure, so the slightly deeper discon- 
tinuity of K8 is probably not significant. The Fuki.o (1977) model ARC-TR has a 
pronounced "bump" in the velocity-depth profile at 500 km because of a rapid 
decrease in p(A) at 20°. GCA also has a very slight inflection near 540 km to 
match strengthened first arrival amplitudes near 21°. While all the models are 
very similar in the 650 km depth ra.ige, GCA is derived from the most complete 
data set in toe 2(>® - 30° interval. The seismograms shown in Figures 4. 5, 6 are 
very sensitive to the fine structure of the 660 km discontinuity as is demon- 
strated in Figure 18. The T? model contains the strong gradient just above 670 
km ’ ut a lack oT convincing data near 27° precludes delineation of the transition 
shape just below the break. Grand ic Helmberger (1983) include increased gra- 
dients from 860 km to 750 km for both shield and tectonic S-wave models, while 
a P-wave model for a shield has a similar form (Given & Helmberger, 1980). 
Dziewonski & Anderson's (1981) model PREM also features an increased gradient 
near 700 km on a global scale. 

The differences between the four models with depth are further illustrated 
by comparison of cumulative one-way vertical travel times. In Figure 21. vertical 
travel times are computed for each model starting at 620 km in 20 km steps and 
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are compared thp Herrin (1968) velocity model. It is evident that the curves 
for T7 and GCA are very similar from 620 km to 300 km, then T7 flattens out with 
respect to the Herrin times while GCA continues to be very slow relative to the 
1968 model. ARC-TR and K8 also have nearly constant values for for 

100 - 200 km depth. The arc and shield models’ residuais decrease above 150 
km as expected from their faster uppermost mantle velocities, while T7 and GCA 
become more positive for that depth range. The various crustal thicknesses 
have large effects on the overall vertical travel times. 

Conclusions 

Earthquakes in Mexico recorded at SCARLET form a unique, high-density 
digital short-period P-wave data set for upper mantle study. We have analysed 
over 1400 seismograms and utilized travel time, apparent phase velocity and 
relative amplitude information to produce a tightly constrained, detailed model 
for depths to 1000 km beneath an active oceanic ridge region, the Gulf of Cali- 
fornia. Unusually low velocities to depths of 350 km characterize the spreading 
center model, consistent with teleseismic P and PP residuals of Rowlett 6c For- 
syth (1979) and Dorbath 6c Dorbath (1981), respectively. The abrupt cutoff in 
distance of observations of the travel time branch (AB) for which rays bottom 
just above 4C0 km leads to an interpretation of an anomalously high velocity gra- 
dient from 225 to 390 km, rather than a large velocity step at the base of the low 
velocity zone. Grand 6c Helmberger (1983) make a similar observation for shear 
waves from East Pacific Rise earthquakes. Thus the ridge model, GCA, differs 
significantly from models observed for shield, young continental, and arc 
regimes to 350 km depth. 

Abundant data from 19° to 28* constrain the detailed shape of the 860 km 
discontinuity. A small, 2.8% jump in compressional velocity is accompanied by 
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incraased gradient? both above and below the break. 

Present available data are not sufficient to confirm or reject the idea cf 
undulations of the major velocity discontinuities in the upper mantle. Lateral 
velocity variations between shields, arcs, and tectonic-continental regions seem 
well established to depths of 200 km. This study demonstrates that oceanic 
spreading centers have velocities slower than the other regions to the even 
greater depth of 350 km. 
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Ffgure Captions 

(1) Location map for this study. Stars are epicenters of the 22 earthquakes 
closer than 30°. Small (2°) po-tions of the great circles between SCARLET 
end the events are also shown, indicating the area covered by model GCA. 
Note that all the arc segments fall within the region affected by the Gulf’s 
spreading. Circles locate the calibration events which are farther than 30° 
from Pasadena. 

(2) Stations of the southern Califonia array u:;ed in this paper. Empirical 
corrections for these 98 stations are shown in Figure 9. Latitude is in 
degrees north, longitude in degrees wesL 

(3) An example of an event record section recorded at SCARLET along with the 
synthetic section predicted by model GCA Distances are not corrected for 
event depth. Amplitudes are scaled to the maximum of each trace. Empiri- 
cal station corrections (Figure 9) have been applied, and the data have been 
filtered with a bandpass of .01 to 5 Hz. For clarity, only a few representative 
seismograms are shown, a) A shallow event on the Rivera Fracture Zone 
(no. 2 in Table l) which shows a weak first arrival followed by the reflection 
from the 390 km discontinuity. t>) Synthetic section for the same event. 

(4) a) A 98 km deep event inland from the Middle America Trench, no. 14 in 
Tabic 1. Note the increasing strength of the first arrival near 20° and the 
strong arrival from the 660 km discontinuity. Also the "back branch” of the 
390 km triplication (see Figure 10b) is not in evidence. See Figure 3 for for- 
mat explanation, b) Synthetic section for event 14 using GCA. c) Synthetic 
section for the same event for model 1*7 (Burdick A Helmberger, 1978). 
Note the strong AB branch extending to 22°. which is not visible in the data. 
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(5) a) Format is the same as Figure 3. No. 17 of Table 1 is 56 km deep. Again a 
strong first arrival and reflection from the 660 km velocity discontinuity are 
visible with no sign of the AB trave 1 time branch past 20° . b) GCA synthetic 
section for event 17. c) T7 synthetic section for event 17. Again the AB 
branch is a very strong phase in these synthetic seismograms, but not 
observed in the data. 

(8) a) No. 27 (Table 1). in the same format as Figure 3. Here the "back branch” 
of the 660 km triplication moves out with increasing distance, and the rela- 
tive amplitudes change from a weak first arrival near 26° to a simple pulse 
near 28°. The coherent energy at 109 seconds is pP of this BO km deep 
earthquake, b) GCA synthetic section for event 27. 

(7) a) Data record section of 10 events spanning 9° to 40°. Amplitudes and 
filtering are as in Figure 3. Station and depth corrections have been 
applied. The mantle triplication phtses are clearly visible, b) Synthetic 
record section for GCA on the same scale, computed for surface focus. 
Source wavelets vary with distance. 

(8) Comparison of event 20 (Table l) before and after application of empirical 
station corrections. Record sections are set up as in Figure 3. a) 
uncorrected, b) corrected. Note the improved alignment of the traces 
near 30.3° .31. 9° .33.1*. and 34.4*. 

(9) A comparison between the teleseismic P residuals of Raikes (i960) and the 
empirical station corrections used in this study. Raikes* values (b) are 
referenced to station CSC and include corrections for elevation, sediments 
and crustal thickness which are not included in (a). Still many similarities 
exist, such as relatively negative values near (35°, -116°). positive residuals 
at (34. 5*. -116.5* ), negative (fast) corrections in the Santa Barbara Channel 
(34*,-120*) and negative values southeast of the Imperial Valley (33* , 
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-115*). Differences occur due to elevation, sediment cover and crustal 
thickness in the Peninsular Ranges (33* ,-116°) and the Ventura Basin 
(34.5*, -119°). The -contour interval is 0.2 s. 

(10) Trave’ time data Tor the 29 events, a) Travel times vs distance . reduced by 
10 km/s and corrected for elevation, depth and receiver structure. De.ta 
coverage is almost continuous from 9° -40°. There are 1753 data points; 438 
are secondary arrivals, b) Models GCA and GCA’ superimposed on the 
baseline-shifted travel times. GCA’ is constructed for the axial region of the 
Gulf of California from only very close (9 # -13 # ) data, while GCA represents 
adjacent continental areas; it has a thin lid and small low-velocity zone. The 
two models merge -below 150 km. Data near (22*, 60 s) correspond to the 
anomalous p - A point of Figure 11. Letters refer to travel time branches 
discussed in the text. 

(11) The 58 p - A points plotted with GCA. Triangles indicate measurements 
made with first arrivals; circles are secondary phases. The data point at 
12.0.s/* t 23° is from event 26 and is discussed in the text. 

(12) The results of tau inversion of the travel times shown in Figure 10. The solid 
line is the inversion of the best tau values, while the dotted lines represent 
error bounds implied by error estimates for each tau point. The single 
in ,- ersion model is used as the starting model for forward computations or 
synthetic seismograms. 

(13) Model GCA. Valid for the Gulf of California spreading region, GCA features a 
20 km crust and low velocities to 350 km depth, with an unusually large 
velocity gradient from 225 to 390 km. Velocity discontinuities are 4.92 at 


390 km and 2.3% at 680 km. 
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(14) The top 150 km or GCA compared with GCA’, the model constructed by tau 
inversion of arrival times from earthquakes closer than 13*. GCA' has no 
low velocity zone and a transition zone at the crust-mantle boundary. The 
difference between these models may represent lateral crustal variations 
between the Gulf of California itself (GCA') and the adjacent continental 
areas. 

(15) A comparison of waveform data for distances of 14* - 18* with various 
models. The top row is the data, followed by GCA. the tau starting model 
and T7 (Burdick & Helmberger, 1978). Distances are corrected for depth. 
All seismograms are from event 2 (Table 1) and are scaled and filtered as in 
Figure 3. At A = 14.7* and 15. 2°. the tau and T7 models predict a first 
arrival much stronger than the 390 km reflection., while GCA comes closer 
to the true relative amplitudes. The models are very similar at 16.3* and 
17.8V 

(18) Record section fcr event 26 (Table 1), reduced by 11 km/s. The line indi- 
cates the arrivals picksd as the AB branch in Figure 7 and used in the tau 
inversion. Notice the poor signal coherency and signal-to-noise ratio com- 
pared to the sections in Figures 3, 4. 5, and 6. 

(17) Same as Figure 15 for 20* to 23*. All data are from event 14 (Table 1) 
except for the trace at 21.7* which is from event 17. At 20.3*, GCA 
correctly predicts times and amplitudes of the two arrivals. The tau model 
and T7 both have large intermediate arrivals which are from rays turning 
just above the 390 km discontinuity. The three arrivals in the T7 synthetic 
are clearly not in the data. At 20.9*. the AB branch arrival is interfering 
with the reflection from the 660 km velocity jump in T7 and the tau model, 
causing the relative amplitudes and timing to be off. All three models 
predict the data at 21.9* rather poorly , but again for the tau and T7 models 



the large phase is that of the AB branch. At 22.3' and 23.1', GCA provides 
the best fit in terms of subtle timing and amplitude comparisons. Here the 
AB arrivals would be more than B s behind the first break. 

(18) Detailed analysis of the 680 km discontinuity. The synthetics are calculated 
for model GCA with differing gradients above and below 660 km. from a sim- 
ple step (tc*p row) to large gradients both above and below (bottom row). 
The relative amplitudes and timing of arrivals are much better for the bot- 
tom model than for the other two attempts. Distances are corrected for 
surface focus. 

(19) Four models for differing tectonic regimes. T7 (Burdick A Helmberger, 
1978) is valid for a tectonically active continental region. Island arcs are 
represented by ARC-TR (Fukao, 1977) and shields by K8 (Given A Helm- 
berger, 1980). GCA (this study) represents an oceanic spreading center. 
KB, T7 and ARC-TR are very similar below 200 km depth, but GCA is substan- 
tially slower than the other models to depth or 350 km. 

(20) A comparison of TNA. Grand and Helmberger's (19B3) shear-wave model for 
Mexico and the western United States, and GCA. Note the great similarity in 
general character between the models. TNA has a 4.73C velocity jump at 405 
km and a 7.8T change at 660 km. 

(21) For a starting depth of 620 km. cumulative vertical one-way travel times 
are computed for four models and compared to the 1968 Herrin mantle 
velocities Symbols are plotted at 20 km intervals. T7 and GCA are very 
close from 600 km to 320 km, then T7's residuals flatten out but GCA's 
become increasingly positive. This illustrates the very slow GCA velocities 


above 350 km. 
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Tables 

(1) The epicentral data are taken from the PDE Monthly Listings of the U. S. 
Geological Survey. For the dT/dA data, F denoies first arrivals, while L 
means later-arriving phases. The dT/dA determination for event 10 is 
unreliable and was discarded. The listed error estimate is that of the least- 
squares plane fit to the travel times. Arrival data are used only for the 96 
stations which have empirical receiver corrections. 

(2) Velocity model GCA. 

(3) Velocity model GCA'. This model, which merges with GCA below 150 km, is 
valid along the axis of the Gulf of California. 



Table 1 


ORIGINAL FAC r.:. 
OF POOR QUALlT 


IYKKT DATA 

Ipleanlial Inbarmatlaa 67/ 46 Data 


■rani 


Data 


Origin Tima 

Latltuda 

Leogltuda 

Dap lb 

MagnKuda 

4 

r 

RMS 

Arrival 

Numbar of 

Rumba r 


Mo 

Yaar 

Hr 

Mia. 

lac 

digram 

dagreel 

kn 

"3 

da gran 

./• 

Error, a 

T>pa 

BUUnna 

1 

99 

r«b. 

197* 

90 

jO 

90 J 

17 36.40 

-101 0.30 

34 

3J 

91.89 

10.76 

.106 

F 

9 












36 Be 

10 33 

JOS 

F 

11 












■1.93 

9.39 

.IK 

L 

■ 1 

■ 

93 

** 

1977 

13 

93 

9 

90 01.93 

-109 964 

36 

S7 

16.75 

USE 

.330 

F 

35 












1663 

11.13 

.191 

L 

■8 

1 

10 

Die. 

1979 

oe 

39 

351. 

>649 

-109 3636 

13. 

63 

9.79 

14.44 

J79 

F 

16 












11.83 

1636 

.315 

F 

14 

4 

99 

Oct. 

1971 

u 

36 

97.1 

•7 36.36 

•loeeore 

SO. 

U 

17.93 

19.60 

J40 

F 

99 












19.03 

10.30 

.137 

F 

11 

» 

99 

Jan. 

1979 

ID 

04 

M.O 

17 94.73 

-100 3669 

41. 

68 

91.36 

10.66 

JB82 

F 

33 












■613 

LO. 06 

J64 

F 

18 

• 

99 

Hen. 

1979 

30 

49 

49.3 

16 11.10 

-96 3760 

SB. 

67 

84.40 

9.92 

JM 

F 

■3 












■(.S3 

611 

064 

F 

■6 












■634 

10 09 

.170 

L 

49 

T 

99 

Dae. 

1971 

93 

97 

(3D 

10 91.96 

-103 3160 

to 

63 

■9 33 

9.00 

.379 

F 

33 

« 

. n 

br 

1979 

01 

39 

14 JO 

17 01.36 

-99 44.10 

36 

63 

■1.94 

10.44 

.160 

F 

n 












■3.44 

693 

J93 

F 

■6 












■9.40 

9 14 

J77 

L 

39 

• 

m 

la 

1971 

OB 

JO 

94.3 

ir a oo 

-94 3064 

107. 

AJ 

■3.33 

9.19 

J099 

F 

41 












■7J3 

9.03 

.064 

F 

91 

to 

90 

M»J 

1979 

11 

B 

41.0 

94 46.49 

-100 03.36 

10. 

66 






It 

10 

Sap. 

1979 

£3 

94 

15.6 

14 16.90 

-91 9969 

•4. 

66 

SO 86 

690 

J031 

F 

08 

1C 

11 

Ian 

1979 

CO 

06 

04 2 

16 46.B4 

-107 06.40 

36 

64 

19.77 

19.61 

.m 

F 

14 












1673 

1LM 

■7 

F 

■ 

19 

19 

Dae 

1979 

or 

39 

19.3 

19 36 49 

-106 36.36 

36 

69 

13 96 

1646 

220 

F 

90 












13.63 

11.39 

265 

L 

10 

14 

99 

Sap 

>979 

te 

91 

41.0 

16 36 90 

-109 16.78 

96 

63 

90.86 

10.73 

■ 1 

F 

17 












■.34 

10.34 

.149 

F 

16 












■0.79 

9.46 

.339 

L 

87 

U 

99 

Fab. 

1979 

09 

19 

17.0 

19 3630 

-100 1606 

86 

6J 

■ 10 

10.90 

.186 

F 

11 












91.81 

10.03 

.107 

F 

11 

19 

09 

Jan 

1979 

u 

31 

30.7 

16 16 BO 

-109 47.10 

36 

69 

19.69 

10.97 

339 

F 

■ 1 












■1.73 

10.94 

2*7 

F 

10 












■SB 

9.43 

>493 

L 

14 

17 

09 

Jul 

1979 

90 

13 

■3.3 

16 99.99 

-100 0 49 

09. 

66 

11.98 

10.46 

.in 

F 

14 












■ 91 

10.09 

.173 

F 

14 












■1.90 

• 63 

J914 

L 

48 

19 

93 

kat 

1979 

DO 

36 

31.2 

10 1694 

K 1630 

96 

67 

97.91 

649 

D46 

F 

77 

19 

n 

Oat 

1979 

14 

33 

97.3 

13 4696 

-90 6966 

•6 

67 

91.44 

633 

J039 

F 

70 

00 

09 

Dae 

1979 

11 

33 

34.0 

13 0670 

-69 3610 

36 

66 

■ 80 

679 

.036 

F 

77 

■1 

97 

Dal 

1979 

91 

43 

94.9 

13 4636 

-90 4630 

06 

66 

91.87 

663 

■039 

F 

74 

99 

91 

M«J 

1979 

01 

07 

914 

19 46.09 

-67 09 4 9 

76 

64 

34 89 

666 

035 

F 

69 

93 

09 

Dae. 

1979 

14 

98 

39.3 

13 49 73 

-09 17.98 

39 

67 

■ 44 

636 

.040 

F 

71 

94 

14 

Mar 

1979 

It 

01 

94.3 

17 87.19 

-101 16.96 

39 

16 

■ 60 

10.36 

.163 

F 

14 












■ 43 

1040 

.187 

F 

■I 












11 tt 

9.44 

.332 

L 

60 

99 

19 

Mar 

1979 

to 

It 

11. T 

17 39 73 

-100 30 46 

36 

64 

■ 60 

10.79 

.133 

F 

■0 












■ 43 

10 49 

183 

F 

17 







Iplotntrol Inform* tloo 


■rant 


Dal* 


OrlflnTImo 

Latitude 

UmfUud* 

Number 

d»t 

No 

Y.»r 

Hr. Min Boo. 

dl|DM 

dogreee 

X 

K 

Oct. 

lire 

14 07 06 

IB 10.96 

-1 CM *8 88 


tr 

4 


1«T» 

08 

■8 

41.7 

18 41.04 

-33 38 78 

m 

10 

Jon 

1879 

IS 

84 

14 J 

1IHJ4 

-83 38.58 

m 

K 

J*a 

1879 

IT 

10 

44.1 

17 39.94 

-100 58.88 


<n/dA Dau 


Depth 

Magnitude 

A 

V 

RMS 

Arrteed 

Number of 

km 

"l 

do (ran 

./• 

Error. • 

Typ* 

Station* 



•1 63 

■ 49 

■49 

L 

38 

39 

U 

81 86 

10.81 

644 

r 

85 



8918 

1003 

■68 

r 

10 



■604 

648 

■ 10 

L 

34 



*1.36 

1608 

■38 

L 

85 

80 

67 

8683 

600 

JD63 

r 

61 



■7 08 

681 

■35 

i 

18 

186 

60 

■601 

600 

J064 

r 

48 

36 

64 

■C 98 

LOTI 

.158 

r 

18 



•687 

10 44 

.147 

f 

18 



■LTB 

656 

J8T 

i 

■6 


on&VAL r, t . M 

OF POOS QUALITY 



OKiiiMAi. r'i r_ 

OF POOR QUALI1Y 


Table 2 



Model GCA 


Depth 

Velocity 

Depth 

Velocity 

km 

km/ 8 

km 

km/s 

0 . 

6.400 

300. 

0.403 

19. 

6.400 

325. 

0.520 

20. 

7.000 

350. 

0.630 

35. 

7.750 

375. 

0.750 

50. 

7.700 

390. 

0.019 

75. 

7.050 

391. 

9.250 

100. 

7.900 

450. 

9.476 

125. 

7.938 

530. 

9.000 

150. 

7.975 

620. 

10.060 

175. 

0.013 

660. 

10.360 

200. 

0.050 

661. 

10.05'J 

225. 

0.100 

600. 

10.760 

250. 

0.160 

970. 

11.340 

275. 

0.285 




i 



Table 3 


n^‘^ AL PAe "- « 

V POOR (&. Model GCA' 


Depth 

Velocity 

ion 

km/s 

0.0 

6.200 

5.0 

6.400 

12.0 

6.700 

25.5 

7.000 

31.0 

7.500 

51.0 

7.750 

100.0 

7.870 

115.0 

7.970 

150.0 

8.000 









ORIGINAL PAGE 13 
OF POOR QUALITY 





ORIGINAL PAGE IS 
OF POOR QUALITY 


Time-A/9, s 

24 34 44 24 34 44 

. , — . I "7 





ORIGINAL fr :C 
OF POOR QUALITY 


Time-A/IO, s 



















OR!G ! !VAL PAG 7.' [J 

OF POOR QUALITY 


82.6 92.6 




102.6 

T 


(q) 12/6/78 Uncorrected 


Time-A/12, s 

112.6 122.6 82.6 92.6 102.6 112.6 122 .6 

fnlyV*»ifWi/yrA*Yty*W\yV\fAW l 


^\VvV\\VA\MA^- 

-^\i\ 

. V .. . «. \ \ » .V «v *»A »V -rt 


31 


A, 


'-) 32 


■A y > 

; I W 

-it j' ' 3 3 

- Vy\i 

-v vv^^'AA , *^V^\' v v\/VV v ' 


_a ^A'^V^/A^V — v-W^WWV- 

— A* 

— . — yV 


i. 

tVi 


\\ a AAV\^~v v ' — 

vv ' , ir 


34 




(b) Corrected 





— T| AjS/\(VV'^~ v ^VL'' N ^/ r VV^^^ 







l' oV v^v^^ v ^^W\A^''n/wv^^a 

I <v^ 

^I^VvVVV^W \W v ^Ar^Vv\VVV^ 




ORIGINAL PAGE 13 

OF FOOR QUALITY 

















Depth, km 


ORIGINAL PAGE !J 

OF POOR QUALITY 


too 


DATA 
Vel.. km/s 
8 9 


300h 


503 L 


100 


jOOJ- 
500 L 


IOC 


300 - 
50qL 


6 9 ^ 


A ■ 14 . 7 ° 
CCA 

w W' A/\/V- 

tou 


T7 

wv — 


15.2° 


16.3° 


— styVWv — 

SYNTHETICS 
— v\/L- -aa/W\A 


— V\A^-^\a- -a yvy^\- 


'5 s' 


mAv\/\ — 


17.6° 

—A/^VW 

-aa/Ua 

-A/\/fy\/\r — 


Fig 15 




ORIGINAL PAGE IS 
OF POOR QUALITY 


Time-A/ll, s 

65 75 85 95 





original nr: hi 

OF POOR QUALITY 



/%. 17 







ORIGINAL . C AG E is 

OF POOR QUALITY 


Velocity, km/s 






Depth, 


ORIGINAL F AG- r S 

OF POOR QUALITY 


Velocity, km/s 



Fig. 20 



ORIGINAL PACE t'3 

0F p °OR QUALITY 


CUMULATIVE ONE-HAY TRAVEL TIMES 


□ 

<n 


GCA 


RRC-TR 




